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SUMMARY 

— tSS ^ i l“ s ssnrr? 1 . . Tr ca i 

^t'; 1 data ” IU be c ”" pa " d ” irt ^.AuLi^suLs 12 iu2^5£; t 

light gliders. Pr ° gram WhlCh may b6 USGd 10 the design and evaluation of ultra- 
begun^n^ry' 3 ^ 8 ;^ 1568 ^ 8 ° alS ^ ° f the investigation 


INTRODUCTION 

In recent years the performance and variety of hang glider designs have 
increased dramatically. Flight conditions and demands that are placed o^ang 
f!“ " ”! Ver l ' d,ft -ft fro. those encountered by older Z H l s Llr'.t * 

elide rI«S ol ITosl to rn'IndT” T ^J 8 ag °’ ^ deslgns achieve 

Srb a 2 2«w h ™i l8b , a ? 60 °? ” <19 'r “•> <Ref - »' -ddiuon n 

bulent) thermal .lying, increased controllability has made limited aerohar-ir 
rn.neuv.rs possible. Several years ago the results of NASA w!„d 

o 2h°S 2.2 s?iS ef ; ^’‘V" used to obt.™ !™ till 

l 1 1 f v ew desl 3 ns • Although not ail flight regimes and 

useful The a hair»lir e 5 horou *J 1 y investigated, the data that did exist proved 
seful. The hang glider has evolved, however, to the point that these original 

bSTslMwHlte" 8) virh' er ! " , ' >lled ; Th « f1 1 E ht characteristics of nodem 
5 to 7 6 22 !2lls 1, 'Im CRtcndlng to 31 n, (36 ft.), aspect ratios from 

data for the high' 

?:r; f d ^s£,r st ' 

such :s ?raX p :r:r,r 

although recovery is theoretically possible in less than 15m (50 ft! R^f Z 
Statistics from hang gliding accidents in 1977 and 1978 show that, despite a moie 
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. v.w t-hp industry in the last few years, such 

SSS.SS’S.’SrS." 2-“ » - — 

Work was begun in January ^f^esi gn and 8 e^lultiorof C modern d hang P gliders, 
quantitative tools for use in the f t nd c i 0S ely integrated phases: 

The investigation consists of ^ w ° ts wil i be made on scale models in 

1) Basic force and moment measurem ^kSPi Ames Research Center. Models 
one of the 2 mx 3 m wind tunnel^ at geometric, elastic, and 

aerodynamic°properties o£ * ruprnsentarlve class c£ -ode™ glider. 

2) A computer program, based ^^e-elLen^strictu^a^methods, for 

potential aerodynamics and aeroelastic corrections 

predicting the ne.scredarondsvth ^ ^ tests> 

Sira" d ;e promulgated for t/e analpsis oi iutnre glider 

designs. 

As this research is to be conducted 
describes the goals and general approach of the prc.ect 
published at a later date. 

WIND TUNNEL TESTS 


Models 

_jsn grfnp“r:; 5 -:«ie c :rdiu « 

full scale value. 

Although there exist today a «U. X'SLSl U “ 

no longer possible to test a ^tanda aircraft, sufficient similarity does 

predict universal characterist c determined from tests on a limited 

exist so that certain characteristics may be determin^ feature8> In 

number of models and applied to many °^ rties g of such gliders ma y be obtained 
this way, good approximations P but carefully selected, geometries, 

from tests on a small group with di * types, from the older Rogallo- 

The models selected span a wi e ra g „ d hig h performance designs. 

ss e«r.n 

llV,lZ[Z S \rZ\Z, on rhis group oi nodcls. 

f pinotic scaling has been demonstrated recently (Ref. 11) • 
The importance of elastic seal g fc{j vary consld erably with 

The flight characteristics of gli flexibility of the frame and defor- 

changes in loading. This is caused by this reason, it is infor- 
mation of the sail of these Jl* „ to remain geometrically 

tant that scale models be constructed in such a way a 
similar to full size gliders under corresponding loads. 

Another key assumption underlying the Rg^This is^because rather complex 
attainment of full-scale Reynolds’ number. Re. This 
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separated-f low effects are anticipated at the larger values of a and B . 
Since available wind tunnels operate at essentially sea-level conditions, it 
follows that any resultant force, F , experienced by the model must equal 
the corresponding F at full scale? (Mach number effects are negligible at 
these "raicrosonic speeds.") 

Force equality can be reasoned from the fact that: 

Re = Re c 
m f 

where the product of speed and typical length, V£, must be the same at both 
scales. With air densities: 


and forces proportional to pV 2 £ 2 , then: 


The combination of equal force and equal strain requirements lead to dif- 
ficulties in the construction of elastically scaled models. Consider that both 
the model and full scale gliders are constructed of tubes and cables of approxi- 
mately circular cross sections of radius r , supporting the fabric sails. 

Since r should be proportional to i for aerodynamic similarity, the strains in 
these tubes are proportional to Fr£/EI or to Fr 2 /EI, with El the familiar 
bending rigidity. The severe requirement on model construction is to ensure 




For a typical glider, assembled of thin-walled aluminum tubes, this quantity is 
of the order 10 -/ N -1 (4 x 10“ 7 lb.* 1 ). If the same construction and 

material were employed on the model, one would get 

/ \ i9.c\ 2 / 2 \ 


(£\ m 2 (£) 

\ ei L \*J \ e 4 


This factor of 25 at one-fifth scale is quite unacceptable. Since weight is 
not believed to be a very significant factor, the situation can be alleviated 
by going to solid-section cylinders of stiffer material on the model. Models 
constructed of steel in this manner approach the desired stiffness; 



It appears that the requirement of equal strains, therefore, can only be 
met by some relaxation of the Reynolds number requirement. The following values 
correspond to the model construction above: 
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This difference is not large and can be reduced further with the use of tubes 
and cables of slightly larger than scaled radii. 

Especially for newer hang glider designs with low billow, it is important 
to duplicate the stretching of the fabric sail as well as the bending of frame 
elements. This requirement may be seen approximately as follows. 

Requiring equal strains in the model and full-scale glider sails, 
for geometric similarity implies that 

o_ _ p (x) dx _ p(x)dx _ F 
E EdA Etdx Eti 

is the same at both scales. Now since j and from above we have let 

F = 1/3.4 F we require that (Et) = 1.8 (Et) . This can be achieved with 
m t m f 

the appropriate choice of Dacron fabric. Values of Et for Dacron sails are 
given in reference 24, from which it can be seen that the proper (Et) may be 

achieved with two layers of material slightly lighter than that used m on 
full-scale gliders. 


Data Reduction 

From measurements of the basic forces and moments on these models, the 
following performance coefficients and static stability derivatives can be 
calculated. 
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Data will be obtained generally at angles of attack, a from -45° through 
stall and at sideslip angles 8 to ± 20°. Tests will be conducted at 
various pressures to obtain data on the presumably significant variation of 
these quantities with the dynamic pressure, q (elastic effects). Test results 
will be corrected for jet blockage and wall effects. 

Much of this data could be used immediately for design purposes with little 
intermediate manipulation. With the use of data on pitching moment coefficient, 
the longitudinal equations of motion may be numerically integrated to show the 
effectiveness of "weight-shif t M control, including required bar pressures (stick 
forces), under various flight conditions. Stall, dive recovery, and other 
aspects of longitudinal motion will be analyzed. A similar analysis for lateral 
motion, taking account of the unusually large coupling between longitudinal 
and lateral modes associated with hang gliders, will also be carried out. At 
the present time, the first wind tunnel model is being constructed at the 
machine shop facilities at Stanford. The frame will have two possible nose 
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angles and by attaching different sails many configurations may be tested. 
Although a final list of configurations to be tested has not yet been deter- 
mined, a tentative group of test models is described in table 1. 


THEORETICAL ANALYSIS 


In conjunction with the testing program, theoretical aerodynamic and 
aeroelastic methodology is being applied toward the development of a computer 
program which will undertake to predict some or all of the quantities measured 
in the experimental portion of the project. 


Several theoretical treatments of "parawing" aerodynamics were published 
in the 1960's (e.g. Ref. 12-14). Lifting surface theory was used to predict 
lift and moment of various parawing configurations with the assumption of a 
particular mode shape (generally taken to be a portion of a right circular 
cone). Induced and profile drags and the effects of rigid leading edges were 
treated. Recent experimental work (Ref. 11) has shown, however, that changes 
in sail shape with angle of attack and dynamic pressure are extremely important 
especially for current hang glider designs. Thus, not only is the assumption 
of conical canopy shape no longer valid, but no rigid analytic assumption of 
mode shape can be used. 




The approach taken in the present analysis consists of two major parts: 

1) The determination of airloads for a prescribed mode shape, and 

2) The flexible structural response to this calculated loading, 
resulting in a new approximation for canopy shape. 

The iterated procedure, shown schematically in Fig. 1, is used to obtain a 
solution for pressure distribution without the need for specifying the exact 
sail shape initially. 


From these predicted airloads, force and moment coefficients may be 
calculated and compared with experimental results. 


Aerodynamics 


Linearized, steady, lifting-surface theory for incompressible flow is used 
in the prediction of aerodynamic loads on the glider. Under such conditions 
the flow over the glider satisfies Laplace's equation: V ^cf> = 0 which 

may be solved with the use of vortex-lattice or kernel-function methods. The 
approach taken here utilizes the former method described by Woodward and 
Rubbert (Refs. 15,16) with a code by Nathman (Ref. 17) used at Stanford's 
computing facilities. 

The sail is divided into finite elements as shown in Fig. 2. Each element 
is idealized as a flat panel of constant doublet strength, K, defined as the 
discontinuity in potential between the upper and lower surfaces, 


K 



t 


* 
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As shown in the appendix, this leads to the following expression for the 
velocity induced at points outside this surface 

{V) - ± [C]00 

where [C] is the aerodynamic influence matrix described in the appendix. The 
doublet strength for each panel is chosen so that the flow at the surface of 
the glider is tangent (zero sail porosity). This condition is satisfied if the 
normal velocity induced by the system of doublets just cancels the free-stream 
normal velocity: 

[C ]{K} = - 4ir {n • V } 

n oo 


Since the surface normals and influence matrix may be computed from the assumed 
sail geometry and since the free-stream velocity is given, the value of K can 
be calculated over the surface. 

Once K is known, the vorticity on the surface is given by: 


*Y = n x V K 

and the loading: 

£p(x) = p V x y 

These pressures are then used to calculate the desired force and moment coeffi- 
cients according to standard definitions (cf. Ref. 18). The procedure is 
summarized in Fig. 3. 

Figs. 4-6 show the preliminary results of this theory applied to some 
simple planforms for which experimental data is available. (Ref. 19,20). Agree- 
ment is close although effects of leading edges and deviation from conical 
geometry are not considered. 


It should be noted that these results are the predictions of the aero- 
dynamic portion of the program only. A rigid mode shape is assumed and so agree' 
ment with experiment can only be expected at intermediate a . The combination 
of this portion of the program and the structural analysis described below is 
presently underway and results are not yet available. 

This analysis does not include the effect of pilot, cable or frame inter- 
ference. It applies only to unseparated flow and does not include viscous 
effects. Corrections to the first stage of the analysis, taking these effects 
into account, are being studied and can, hopefully, be implemented in later 
work. 
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Structural Analysis 

assm^dtThl ?n d °V han f 8llder constitute a ra <*er flexible structure 

loaded t a state of quasi-static equilibrium. Tension members, axially 

clefrlv b dpf?n 5 6nd J n ? members ’ and membrane surfaces are all involved, with 
! fj y r J ! d modal connect i°ns. It is evident that the finite-element 
and k i° j at *-£, structural analysis is the only feasible way of representing 
and balancing the complete system of internal and external loads. 8 

. 1ac P 16 a PP r °ach taken here involves an analysis of the glider frame by 
classics 1 methods and modelling of the sail as a membrane with very small 
flexural rigidity. The procedure is diagrammed in Fig. 7. 

. A a incremental loading technique as described by Turner et al. (Ref 21) 
is used to predict the response of the entire structure to the given applied 
load. The pressure distribution given by the aerodynamic portion of the 
analysis is broken down into small increments and the change in shape due to 
this incremental load !s calculated. This Is done by expressing the pressure, 

fk 1 ’. ,' e3Ch panel ln terms of et i ulva l e nt nodal loads, F., and calculating 

e displacement, D^, of the nodes by tne relation: 1 

{F} = [S] {D} 


^ ^ 3 St ^ fness matrix, made up of a linear, elastic part [S ] 

of tlt^ f f ? r sai1 ^retching (despite the anisotropic stretching behavior 
textile materials, the glider sail is assumed isotropic for the early stage 
° the investigation) and a non-linear geometric part [S„] which depends 
on the geometry and initial tension. The addition of this leometric stiffness 
to the conventional stiffness matrix allows the non-linear strain-displacement 

. atlons associated with this large displacement problem to be incorporated 
in an approximate manner . 

A method described by Argyris (Ref. 22) is adapted here to generate the 

felltin^shln Tv tr± ?' ThiS meth ° d assumes a linear strain-displacement 

techni hlP K ^ hln he elements and is considerably simpler than conventional 
techniques which require calculation of the strain energy (cf. Ref. 23). 

At each step the geometric stiffness matrix is undated and nodal forces 
and incremental displacements calculated. After the step-by-step process is 
C ^ P Kj ed ’i. the inc ; remental displacements are summed to obtain a new mode shape 
w c is then used as input to the aerodynamic program for another iteration. 

A code based on this approach has been developed and is presently being 
checked by comparison with test cases for which analytic solutions are 
possible. Preliminary work indicates agreement to within a few percent in 

displacement although further work is needed to assure convergence in some 
cases • 
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J° r “° n f ° f fj' *— «ig.tlo„ will be used to 
incorporate the effecM of Valldity , md uUl 8““« effort, to 

other phenomena in the analytical potion of”h" C r; s «rci” g ' ed8e 


CONCLUDING REMARKS 


of i; provide * > w <* «. 

the analysis can properlv take account of^l ^ ^ UCh W ° rk is rec ? uired before 
elastic effects aLLt^Uh f ^ 

writing , the aerodynamic and structural rn!tin!!i At ^ of thls 

it is expected that this will ' , ines Klve not been combined although 

presently l« , t , ^ ,^'1" models are " 

from both Che theoretical^ and ewHm.nr ? Uw «<•»• year. Results 

lished as they become available. r> ' 1tt “ ° ' us ri!so »rch will be pub- 


APPENDIX 


Aerodynamic Influence Matrix Calculati 


on 


Expressing the velocity perturbation potential MPl ,► 

toms of the value of * and its normal derivative A on til in , 'i ‘"a 
bv (.reen s theorem: --- ’ on fluid boundary 

^ n 


HP) 


•/“P (6). 


where r is the distance between P and P* . „ • _ , , 

1 r arul 1 i point on the boundary, S. 

order lint" the' f tow iftanJont'''to lhc‘"™faee’' "" r,IHrc ’ then "» * * ”( “" d . *• 


so; 


and 


tangent 

to the 

u 

He 

PhnT = ” 

n 

u 

e 

♦(D = 1 

r 

G.s - 

J 

li 
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If K is assumed constant over each of the panels S., then: 

The aerodynamic influence coefficient of the region S, on the point P is 
thus defined as: 3 1 



Expressing this velocity at several points in matrix notation: 

(V> = ~ [C] (K) 
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LIST OF SYMBOLS 


Notation 


[C] Aerodynamic influence matrix (see appendix) 

D Displacement of panel nodes 

E Elasticity constant 

El Bending rigidity 

F Force 

K Doublet strength 

i Typical length 

n Unit vector normal to surface 

P Point on surface of sail 

P Pressure 

p(x) Loading on sail per unit length 

q Dynamic pressure 

Re Reynolds number 

t Sail thickness 

[S] Stiffness matrix 

V Fluid velocity 

V Free-stream velocity 

a Angle of attack 

3 Angle of side slip 

y Vorticity 

e Strain 

p Fluid density 

a Stress 

(p Velocity perturbation potential 

C T Lift coefficient 


i 


Drag coefficient 

Pitching moment coefficient (based on keel length and referred to the 
c r /2 point) 

Slope of pitching moment curve with respect to a 

Effective dihedral (rolling moment coefficient due to yaw) 

Yawing moment coefficient due to sideslip 

Subscripts 


elastic 
full scale 
geometric 

indicies refer to individual panels 

lower surface 

model 

normal component 


upper surfaces 


TABLE 1 


Details of Proposed Models 


Config- 

uration 

Number 

* 

Airframe 

Basic Design Features 

Comments 

i 

A 

High sweep, low aspect 
ratio, "standard" 

For comparison with more 
recent designs and pre- 
vious wind tunnel studies 

2 

B 

High sweep, medium 
aspect ratio, 2° 
billow "intermediate" 

Comparison with standard 
and high performance 
designs; effects of "billow" 

3 

C 

High performance 
medium sweep (35°) 
zero tip chord 

Washout not fixed by tips 

4 

C 

Same as # 3 with 45° 
sweep 

Effect of sweep on 
stability 

5 

D 

High performance low 
billow fixed minimum 
twist with "floating" 
ribs at tips 

Effect of this common tip 
geometry on C^ 

6 

E 

High performance low 
billow high twist 
dihedral 

Features common to many 
contemporary hang gliders 

7 

E 

Same as #6 with 
decreased twist 

Effect of twist on per- 
formance and stability 

3 

E 

Same as #6 without 
geometric dihedral 

Dihedral effects on lateral 
stability and control 
response 

9 

E 

Same as ^6 with "keel 
pocket" and large reflex 
at root chord 

Reflex effects on longi- 
tudinal stability and 
lateral control 

10 

F 

Similar to #6 with low 
taper planform, low 
twist, reflex 

Common to some of the 
highest performance gliders. 


*Some configurations can be changed with minor model modifications, which 
results In the need for only 6 airframes for the 10 configurations listed. 
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3.- Algorithm for aerodynamic analysis 
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Figure 4,- Results of aerodynamic program 
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Figure 5.- Results of aerodynamic program. 
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Figure 6.- Results of aerodynamic program. 
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Figure 7.- Computational procedure — deflection analysis. 
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